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ABSTRACT 

We propose a novel method to estimate M^/A/bh, the ratio of steUar mass (A/*) to 
black hole mass (AfBH) at various redshifts using two recent observational results: the 
correlation between the bolometric luminosity of active galactic nuclei (AGN) and the 
star formation rate (SFR) in their host galaxies, and the correlation between SFR and 
Af, in star-forming (SF) galaxies. Our analysis is based on Afen and Lboi measure- 
ments in two large samples of type-I AGN at z~l and z~2, and the measurements of 
M^/Mbh in 0.05<z<0.2 red galaxies. We find that A/^/A/bh depends on A^bh at all 
redshifts. At z~2, Af*/AfBH'-280 and -40 for AfBH = 10^ and AfBH=lO^Af0, respec- 
tively. Af*/AfBH grows by a factor of ~4-8 from z~2 to z=0 with extreme cases that 
are as large as 10-20. The evolution is steeper than reported in other studies, proba- 
bly because we treat only AGN in SF hosts. We caution that estimates of A/^/A/bh 
evolution which ignore the dependence of this ratio on AfsH can lead to erroneous 
conclusions. 
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1 INTRODUCTION 

Study of the co-evolution of Active Galactic Nuclei (AGN) 
and their host galaxies provides important clues about the 
growth of super massive black holes (SMBHs) and the 
star formation (SF) history of the Universe. In the local 
universe one finds a tight correlation between the SMBH 
mass (Mbh) and the mass of the bulge of its host, A/buigo 
(Marconi & Hunt, 2003; Haring & Rix 2004, hereafter 
HR04), or alternatively with the stellar velocity dispersion 
a* (Ferrarese & Merrit 2000; Gebhardt et al. 2000). Typi- 
cally Afbuige/AfBH— 500-1000. MBuigc/AfBH must have been 
smaller at high redshift. For example, Mbh ~ io^~^° Mq are 
often observed at z=3-6 (e.g. Netzer 2003; Shemmer et al. 
2004; Fan et al. 2006), yet galaxies that are 500-1000 times 
more massive are never observed at z>0.5 and are very rare 
even at z~0. 

The evolution of the Mbh — o"* relationship has been 
studied in numerous papers. Examples are Shields et al. 
(2003) who find no evolution up to z~2, and Woo et al. 
(2008) who suggest that MBuigo/AfBH has increased by a 
factor of ~3 since z=0.6. The uncertainties in all such mea- 
surements are very large due to the difficulties in measuring 
CT« or Afbuigc in high redshift AGN hosts. While measuring 
Afbuige in high-redshift galaxies is difficult, the total stel- 
lar mass, Mt is easier to obtain. This is achieved by multi- 
band spectral energy distribution (SED) fitting, which is 
used to constrain Af, /L*. However, obtaining Aft/AfBH for 
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AGN hosts is severely limited by the problematic subtrac- 
tion of the bright point-like continuum in type-I AGN and 
the estimation of Mbh in type-II AGN. 

Several recent studies used deep imaging to estimate 
M*/AfBH (or I/Hoat/AfBH) in z~0.5-3 AGN by resolving the 
host galaxy emission and by careful PSF modelling (e.g. 
Kukula et al. 2001; Peng et al. 2006; Kotilainen et al. 2007; 
Bennert et al. 2010). The measured host luminosity was 
translated to i\f, by assuming a certain M^/Lt ratio. A 
common assumption is that AGN hosts are "red and dead" , 
with a stellar population which evolves passively from a 
high formation redshift, e.g. Zform —5. A detailed study of 
this type, including a summary of many earlier findings, is 
given in Decarli et al. (2010; hereafter DIO) who find that 
Af*/AfBH evolves following {M^/Mbh) oc 2"°-^®. The as- 
sumption of non-SF AGN hosts may describe some objects, 
but many studies find that hosts of luminous AGN often 
contain much younger stellar populations (e.g., Kauffmann 
et al. 2003; Jahnke et al. 2004; Silverman et al. 2009; Merloni 
et al. 2010, and references therein). 

In this Letter we suggest a novel method to estimate the 
evolution of M^/Mbh- Our approach makes use of the known 
relationships between star formation rate (SFR) and Af, in 
SF galaxies (SFGs), and between the bolometric luminosity 
of AGN (Lboi) and the SFR in their hosts. We describe these 
relationships in ij2]and use them to estimate M^/Mbh in SF 
galaxies at z~0.15, z~l and z~2 in f|3l In ^we discuss the 
implications to the co-evolution of SMBHs and their hosts. 
Throughout this work we assume a standard ACDM cos- 
mology with nA=0.7, Qm=0.3 and i/o=70kms"^ Mpc~^ 
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GOODS field. E07 studied blue {U - g 
z=0.8-1.2 and found SFR ~ 57 (M,/10 



2 SFR-M, AND SFR-Lboi CORRELATIONS 

Our work is based on two empirical correlations. The first 
is the well established, redshift dependent, SF sequence 
(a correlation between SFR and total stellar mass, M*) 
in SF galaxies. There are numerous papers on this issue 
and the following is only a partial list which is most rele- 
vant to our work. Brinchmann et al. (2004; hereafter B04) 
studied a large sample of low redshift SFGs in the Sloan 
Digital Sky Survey (SDSS). They show a clear relation- 
ship of the form SFR ~ 8.7 (M./10")°'^^ The result 
was later confirmed by Salim et al. (2007) who used a 
combination of SDSS and GALEX data. Similar relation- 
ships at higher redshifts, based on mid-IR and UV ob- 
servations and multi-wavelength SED modelling, are re- 
ported in Elbaz et al. (2007; hereafter E07), Daddi et al. 
(2007; hereafter D07), Noeske et al. (2007), Drory & Al- 
varez (2008) and several other papers (see a more complete 
Ust in Dutton et al. 2009). E07 and D07 used the com- 
bined HST, Spitzer and ground-based photometry of the 

^ 1.45) galaxies at 
iM0)°-'M0yr-i. 
D07 studied feii'-selected galaxies at z~l-3. For z~2 they 
find SFR ~ 200 (M*/10")°'* Mq yr'V The scatter around 
both relations is ~0.3-0.4 dex. Noeske et al. (2007) reports 
SFR ex Af°-^ for z=0.2-0.7 AEGIS galaxies. Drory and Al- 
varez (2008) studied the FORS Deep field up to z~4.5 and 
find similar trends, albeit with systematically lower SFR. 
Much of the differences between the various SFR-M* rela- 
tions can be attributed to the selective inclusion of extremely 
low-SFR galaxies in the samples under study. 

There is mounting evidence that many, perhaps most 
AGN hosts are actively forming stars. This relates to the 
issue of whether such hosts are "blue", "red", or "green 
valley" sources (e.g., Brammer et al. 2009 and references 
therein). This is especially important at high redshift due to 
the known tendency for the fraction of blue galaxies to in- 
crease with redshift (e.g. E07 and references therein). Works 
like Salim et al. (2007) show that low-z AGN hosts popu- 
late the more massive part of the SF sequence. These high 
masses, in turn, mean that the apparent "green" colours of 
many AGN hosts are the consequence of low specific SFR 
(SSFR), not a low SFR. Moreover, Brammer et al. (2009) 
show that many green valley galaxies belong to the blue se- 
quence once reddening is properly taken into account. While 
such ideas are well supported for low redshift AGN, there is 
a need for more evidence at high-z. Some such data already 
exists, e.g. the Silvermann et al. (2009) work that claims 
that the SFR and M. in AGN hosts at z~0.7-l are indistin- 
guishable from those of inactive SF galaxies. Given all the 
above, our first assumption is that the hosts of most lumi- 
nous, optically-selected AGN are part of the SF sequence at 
all redshifts. This assumption is further justified in ^3.11 

Our second assumption is that there exists a signifi- 
cant correlation between the bolometric luminosity of AGN 
(I/boi) and the SFR of their hosts. This correlation is hinted 
at in various intermediate and high redshift studies (Net- 
zer et al. 2007; Lutz et al. 2008) and in low redshift type-II 
AGNs (Netzer 2009; hereafter N09). Here we adopt an up- 
dated version of the correlation presented in N09, by exclud- 
ing LINERs from the group of low-z AGN. This gives 
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Figure 1. The Lbol-SFR relation used in this work. The magenta 
symbols are z~l and z~2 zCOSMOS sources described in the text 
and the small green dots are SDSS S2s. All other subsamples are 
described in N09. The solid line represents Eq. [T] 



SFR ~ 32.8 (Lboi/lO*'^ ergss 
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This equation is not a fit to the data but rather a line that 
goes through the points. As discussed in N09, there is no 
simple way to derive a best fit function to this inhomoge- 
neous data set. 

Since the above correlation is crucial for our study, we 
show the data used here in Fig. [l] In addition to the samples 
presented and discussed in N09, this version also includes 
two additional groups of zCOSMOS type-II AGNs at z~l 
and z~2, kindly provided by V. Mainieri. These are X-ray 
selected sources where Lboi is estimated from L(2 — lOkeV) 
and the SFR is based on multi-wavelength SED fitting. The 
two groups lie on the above relationship with a scatter in 
SFR of about 0.4 dex, comparable to the overall scatter. 
There are several reasons for this scatter. At the low-Lboi 
end, the scatter is related to the inaccuracies in SFR and 
Lboi determination as well as real scatter in these properties. 
At z~l and z~2, the scatter refiects the uncertainty in SED 
modelling and the range of M, across the SF sequence. At 
the high-Lboi end, much of the scatter is due to observational 
uncertainties, the incompleteness of the high redshift AGN 
samples and, perhaps, the extreme SFR in mergers. A more 
complete account of these issues can be found in N09. 

The following analysis relies on the combination of the 
above correlations for SF AGN hosts. Given Lboi, we can 
determine SFR (Eq. [1} and this, in turn, can be translated 
to M* given the redshift dependent SFR-Af* correlations. 



3 THE REDSHIFT EVOLUTION OF M./Mbh 

3.1 M,/Mbh at low and high redshift 

To explore Mt/MsH, we define several samples of both 
AGN and non-AGN galaxies. The first is a large sample of 
0.05<z<0.2 galaxies fi-om the value-added MPA/JHU SDSS 
DR4 database. This includes photometry, M*, SFR and a, 
for all sources (B04). To justify the use of the Mbh-o"* rela- 
tion of Tremaine et al. (2002), we choose only red galaxies, 
following the colour cut of Baldry et al. (2004). The number 
of such sources, after applying some basic quality criteria, is 
210,158. They cover the range 5 x 10^ < M./Mq < 6 x lO" 
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and 10*^ < Mbh/Mq < 2 x 10^ (the latter lower limit was 
chosen to filter out dubious measurements). M* is derived 
from the SDSS photometry and thus the flux limit of the 
sample (rpotro < 17.7) affects the number of low-M* galax- 
ies. Here we deal mostly with the larger M, systems and 
these limitations do not affect our general conclusions. 

The work involves also four AGN samples: type-I and 
type-II samples at 0.1<z<0.2, a type-I sample at z~l and 
a type-I sample at z~2. The first two provide little new in- 
formation but serve to test and to justify the general new 
method presented here. The 0.1<z<0.2 redshift range is cho- 
sen to enable proper measurements of type-I and type-II 
AGN (see details in N09). The type-II AGNs are "strong 
AGN" (Seyfert 2s - S2s but not LINERs) drawn from the 
local SDSS sample. Lboi for these 1152 AGN is determined 
by the [O l] and [O ill] method of N09. The 2814 type-I 
AGN are drawn directly from the SDSS/DR7 (Abazajian 
et al. 2009) database and are analyzed in a way similar to 
the one presented in Netzer & Trakhtenbrot (2007; here- 
after NT07). Here Mbh is estimated from the monochro- 
matic luminosity at 5100A (L5100) and the FWHM of the 
H/3 line. Lboi is determined from I/5100, using the bolometric 
correction factors of Marconi et al. (2004). To avoid selec- 
tion biases, we applied a common observed flux limit for 
the type-I and type-II samples, such that only sources with 
F([Oni])> 4 X 10~^® ergscm"^ s~^ are included. Assuming 
a typical bolometric correction of 2500 (Netzer et al. 2006, 
NT07), this translates to Lboi > 2.5 x lO"^ ergss~\ 

We constructed the SF sequence for our type-II low red- 
shift sample and compared it with the SF sequence presented 
in B04. The agreement is very good, with about 70% of 
the AGN hosts lying on the sequence. We also confirm ear- 
lier results (e.g. Salim et al. 2007) which suggest that AGN 
hosts concentrate at somewhat higher A/* . We find that de- 
spite their apparently red colours, the majority of the red 
S2 hosts are actively forming stars and are situated on the 
SF sequence. These SF S2 hosts are, in fact, mostly green 
valley sources, with lower SSFR. Some 30% of the AGN are 
not part of the SF sequence. We suspect that the fraction of 
such sources at higher redshift is smaller but have no way to 
check it qualitatively. The remaining of the paper and the 
results concerning M, /Mbh refer only to those AGN that 
are on the SF sequence. 

Next, we translate Lboi to SFR for the 0.1<z<0.2 type-I 
AGN using Eq.[T] This is then converted to M, using the low 
redshift SF sequence of B04. Combining with the measured 
Mbh, we obtain M,/Mbh for all these sources. This can be 
compared with M«/Mbh measured directly for the type-II 
z=0. 1-0.2 AGN. The overall good agreement justifies the 
use of a similar procedure for higher redshift AGN samples. 

We selected two higher redshift type-I AGN sam- 
ples from the SDSS/DR7. We chose sources in the range 
0.9<z<l.l and 1.8<z<2 (5330 and 4352 objects, respec- 
tively), which are the nominal redshift ranges of the SFG 
samples of E07 and D07. Full details of the sample selection, 
line fitting and related analysis will be given in a forthcom- 
ing publication. In short, we use a similar procedure to the 
one described in NT07 and in Shen et al. (2008) to mea- 
sure the Mg II A2798 emission hue complex. The FWHM of 
the line and the adjacent continuum luminosity (L3Q00) are 
combined to estimate Mbh using the relation of McLure & 
Dunlop (2004). We calculate Lboi from Z/3000 by calibrating 



iaooo against Z/5100 in a separate sub-sample where the two 
continuum bands are observed in the spectrum. We then es- 
timated the SFR by using Eq.[T] The distribution of inferred 
SFRs (not shown here) clearly shows that these AGN hosts 
are actively forming stars at rates that are comparable to 
the non-AGN SFGs at those redshifts. In particular, 99% of 
the z~2 sources have SFR ^ 42M0 yr~^, the median SFR 
in the SF sequence of D07. 

While most (not all, see earlier comments) high redshift 
AGN hosts are expected to lie on the SF sequence, the range 
in SFR and in M» can be very different from those found 
for non-AGN samples because of the differences in prop- 
erties of the observed samples, in particular different fiux 
limits. To examine this in detail, we focus on our z~2 type- 
I AGN sample. We compare the range of derived SFR and 
possible range of A/, to the same properties in the D07 sam- 
ple, using data kindly provided by Emanuele Daddi. Such 
a comparison involves two crucial factors. First, the D07 
sample includes much fainter sources. As explained above, 
almost all our type-I AGN occupy only the upper part of 
the z~2 SF sequence. This is a direct consequence of the 
SDSS fiux limit. Second, one can consider two approaches 
to deduce typical values of M» for the z~l and z~2 sam- 
ples, by either (1) converting each individual SFR to M* 
through the best fit SF sequence, at the appropriate red- 
shift (i.e. E07 and D07), or (2) sampling the distribution 
of M, per given SFR (in the high-redshift SF sequences), 
for each deduced value of SFR. By definition, the latter will 
result in a considerably broader distribution of the derived 
M« /Mbh , due to the wide range of properties in the observed 
samples. In the following analysis we only use the first ap- 
proach, i.e. we derive M* for each z~l and z~2 AGN host, 
by converting its (derived) SFR through the E07 and D07 
relations, respectively. We note that the ~0.4 dex scatter in 
these relations is a real uncertainty on our high-redshift re- 
sults. The z~2 sample covers 9.5 x lO'^ ^ Mbh/Mq < 5 x lO" 
and 1.7 x lO"' < M,/Mq < 1.8 x lO", while the z~l sam- 
ple covers 3.8 x 10^ < Mbh/Mq < 2.5 x 10^ and 2.8 x 
10^" < M»/Mq < 2.9 X 10". Clearly, the varying ranges of 
Mt and Mbh in the different local and high-redshift samples 
prohibit a simplistic comparison of the mean M«/Mbh. In 
what follows we thus preform a more careful comparison. 

3.2 A'/* /Mbh evolution 

Standard galaxy evolution scenarios suggest that the end 
phase of many high redshift SF galaxies are massive, red 
ellipticals. Therefore, we use our data to compare the prop- 
erties of the z~l and z~2 AGN hosts to those of red 
galaxies in the local Universe. Fig. [2] shows the entire 
sample of 0.05<z<0.2 red galaxies, as a gray scale den- 
sity map. The galaxies form a well-defined band in the 
Af»/MBH-A'/BH plane that follows the approximate relation- 
ship Af»/A'/BHOcMBH~°'' . Also shown is the sample of 
0.1<z<0.2 type-II AGN that follows a similar trend. 

Fig. [2] also shows data for the 30 local galaxies from 
HR04 with dynamically measured Mbh. In this case we 
plot Mbuige rather than M*. All Mbh are taken directly 
from HR04 except those of M 87, NGC 4649 and NGC 4697 
where we used new measurements reported in Gebhardt & 
Thomas (2009), Shen & Gebhardt (2010) and ForesteU et al. 
(2010). This sequence extends up to Mbh ~ 6 x 10® M© and 
Af,/AfBH'^100 and shows a clear dependence on Mbh, sim- 
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Figure 2. Mt/M-^n for different loeal and higfi redshift samples: 
local red galaxies (gray scale density map), type-ll AGN (red 
points), the HR04 sources (green circles) and z~2 type-I AGN 
(blue points). The arrows demonstrate simple scenarios where 
only M* or both M, and A/bh grow by a factor of 2. 



ilar to the red galaxies in our sample and in contradiction 
to to claims of a constant M*/Mbh~500-1000. Several of 
the HR04 galaxies are, in fact, disk+bulge systems with rel- 
atively small bulges and the replacement of Mbuige by M* 
will enlarge the plotted M,/Mbh- Based on this fact and 
the large uncertainties on the estimated Mt/MBH, we sug- 
gest that our sample of local red galaxies is consistent with 
the HR04 sequence. We also show in Fig. [5]M«/Mbh and 
AfBH for our z~2 sample. The high redshift objects form a 
similar sequence, shifted down by about 0.75 dex, with ap- 
proximately the same slope. The z~l sequence (not shown) 
is situated between the two. 

Fig. [2] provides important clues about the redshift evo- 
lution of Mt/MBH- First, A/»/Mbh evolution must be con- 
sidered for a given A/bh. The population mean of A/./A/bh 
would fail to account for the fact that the way such samples 
are drawn is biased towards larger Afen at higher redshifts. 
To demonstrate this, we show in Fig. [3] three histograms, 
representing three vertical cuts in Fig. [2] for 0.2 dex wide 
A^bh bins centred at 10*, 10*'^ and lO" Af©. These cor- 
respond to the case of M^ hut no Mbk growth. The typi- 
cal growth factors of Af*/AfBH between z=2 and z~0.1 are 
~3.8 for AfBH=lO*A^0, ~5.7 for AfBH=lO*-^A^0, and ~7.8 
for A/BH = lO^Af0. We also show, in dashed vertical lines, 
the mean A/./A/bh in the HR04 sample for the same A^bh 
obtained from the best (BCES bisector) linear fit to their 
data. The corresponding M,/Mbk growth factors are ~2, 
~3.3, and ~5.6, for the same values of A/bh. 



4 DISCUSSION 

The analysis presented here suggests strong evolution in 
Af*/A/BH up to z=2, much steeper than what is suggested in 
other studies. In particular, our mean values of A/./A/bh at 
z~l and z~2 are smaller than the ones reported in DIO, 
which is the most up-to-date compilation of such works. 
There are three main reasons for these differences: 

(i) As already mentioned, most earlier works estimate M^ 
by assuming a passive evolution from a high formation red- 



Figure 3. The distributions of A/, /A/bh at different redshift, for 
the three A/bh sub-groups discussed in the text. 



shift. This is obviously an over-simplification of the evolution 
of most galaxies and significant epochs of SF at z<4 (e.g. 
D07; E07; Drory & Alvarez 2008; van Dokkum et al. 2010, 
and references therein) . Accounting for younger stellar pop- 
ulations would result in lower Af« and thus lower A/* /A/bh 
(see DIO and Peng et al. 2006). 

(ii) The DIO sample represents the minority of the AGN 
population, as hinted by two biases. First, the majority of 
the DIO sources lie in the top 15% of the A/bh distributions 
corresponding to their redshift. Second, the selection crite- 
ria for the DIO HST observations could have been biased 
towards large, resolved galaxies with large Af». 

(iii) The mean A/bh in high redshift AGN samples is sys- 
tematically larger than the corresponding low redshift A/bh. 
For example, in our large sample of red galaxies, 99% of 
the sources show AfBH< 10*** A/© while 20% of the z~l 
and 57% of the z~2 AGN have larger A^bh. All the z~2 
sources in DIO have A/bh ^ 10*'^A'/q. Such objects should 
only be compared with local galaxies which host BHs that 
are at least as massive. As Fig. [5] shows, this corresponds to 
A^./A^BH=100-200, instead of the commonly used ~700. 

In conclusion, while our work applies to most AGN, the DIO 
sample probably represents the remaining sources. 

The results presented here point to a scenario where 
many galaxies have to increase their mass by factors of 4-8 
(2-4) since z~2 (z~l). The growth factors for the most mas- 
sive BHs are not well determined since the number of very 
massive galaxies in the local universe is not large enough to 
reliably extend the results of Fig. [2] beyond A^bh — 10^ A^q. 
These numbers represent the requirement for the galaxies to 
ot;er-grow their SMBHs by the above factors. This seems to 
be consistent with models which suggest that the high mass 
SMBHs observed at z~ 2 could have accumulated most of 
their mass by that redshift (Marconi et al. 2004). On the 
other hand, it may be in contradiction with at least some 
scenarios linking AGN activity to the shut-down of SF in 
their host galaxies (see Somerville et al. 2008; Cattaneo et 
al. 2009 and references therein). 

While a full discussion of the various growth scenar- 
ios of A/« is beyond the scope of this Letter, we comment 
briefly on some of these ideas. Major galaxy mergers would 
increase both A/, and A^bh, either through starbursts and 
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gas accretion in "wet mergers" or the possible coalescence 
of the two SMBHs involved in "dry mergers" . However, the- 
oretical and observational arguments (e.g. Lotz et al. 2008; 
Genel et al. 2009) suggest a low rate of such events for z< 2 
galaxies. Thus, major mergers cannot change M»/Mbh by 
more than a factor of ~ 2 — 3 between z=2 and z=0. Small 
"dry mergers" may help. For example, Naab, Johansson & 
Ostriker (2009) show that present-day massive red ellipticals 
gain ~ 40% of their mass through accretion of smaller com- 
panions since z~ 2. Intense SF in outer parts of galaxies due 
to external source of cold gas which does not find its way to 
the centre (e.g. van Dokkum et al. 2010), is another possi- 
bility. More possibilities and more references are discussed 
in Benson & Devereux (2009). The M,/Mbh distributions 
presented here suggest an increase in M,/Mbh by factors 
beyond what is suggested in many theoretical studies. 

Finally, we comment on the possibility that the sug- 
gested evolution of M^/Mbh could be due to two wrong 
assumptions. First, many more AGN hosts may not lie on 
the SF sequence or may not obey the Lboi-SFR correlation 
used here. This is unlikely to be the case at low redshift, 
where SDSS type-II AGNs are used. However, the selection 
of at least some of the most luminous, high redshift high-Lboi 
sources in Fig.[T]may be biased towards high FIR luminosity, 
high SFR AGN hosts (e.g. Zheng et al. 2009) in particular 
if these are found in mergers that are not part of the SF se- 
quence. Herschel observations of well-defined AGN samples 
are likely to resolve this issue. Second, the Drory & Alvarez 
(2008) work shows a decline in SFR at the high-M* end for 
z<2 galaxies. Thus, some of our Lboi-based SFR estimates 
might be associated with considerably larger values of M, . 
The E07 and D07 SF sequences do not show such a decline. 

We conclude that there is a steep evolution in M,/Mbh 
from z~ 2 to z=0 for SF AGN hosts. This trend is barely 
consistent with some, but not all galaxy and BH evolution 
models. We have also demonstrated the crucial importance 
of considering different Mbh groups separately when evalu- 
ating the M*/Mbh evolution. 
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